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Abstract 
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We study the effect of a magnetic field on the interacting quark matter and apply 

o, 

to strange star. We consider the low temperature approximation to strange matter. We 

a ■ 

find that the interacting strange quark matter is more stable compare to free quark gas 
in presence of strong external magnetic field with zero and finite temperature. We then 



calculate strange star structure parameters such as mass and radius and find that the strange 
star is more compact for interacting quark matter than free quark matter in presence of 
strong magnetic field. 
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magnetic field: general 



It has been argued that there is a strong magnetic field in the vicinity of astrophysical objects 
such as neutron stars, white dwarfs and supernovae. Particularly, the magnetic field is very high 



in the surface of the neutron stars. The estimation of the magnetic field strength at the surface 
of neutron stars is done in theoretical models of pulsar emission, (Ruderman 1972) the accretion 
flow in the binary X-ray sources (Ghosh and Lamb 1978) and observation of features in the 
spectra of pulsating X-ray sources e.g., cyclotron lines (Trumper, J. et al. 1978, Wheaton et al. 
1979, Gruber et al. 1980 and Mihara et al. 1990). It has been observed from a sample of more 
than 400 pulsars, the surface magnetic field strength is in the range 2 x 10 10 G < H < 2 x 10 13 G. 

Very recently, two different physical mechanism leading to an amplification of some initial 
magnetic field in a collapsing star have been proposed by several authors (Duncan and Thompson 
1992, Thompson and Duncan 1993 and Bisnovatyi- Kogan 1993). In the first scenario (Bisnovatyi 
-Kogan 1993), the field amplification leads to the formation of an additional toroidal field from 
the poloidal one by twisting of the field lines due to rapid differential rotation of proto- neutron 
stars. The induced toroidal magnetic field could be as large as H ~ 10 15 — 10 17 G, after the 
first twenty seconds of the life of a new-born neutron star. The second mechanism (Duncan and 
Thompson 1992, Thompson and Duncan 1993) is a dynamo action in a differentially rotating 
and convective young neutron star, which is responsible for the strengthenning of some initial 
dipole field upto values of H ~ 10 12 - 3 x 10 13 G or to H ~ 10 14 - 10 15 G. The former value is 
due to the convective episodes arose during the main- sequence and later value is so, if the dipole 
field is generated after the collapse. So, the new born neutron stars have values of magnetic field 
as strong as H ~ 10 14 — 10 16 G, or even more. In the interior of neutron star, it probably reaches 
~ 10 18 G. Therefore, it is advisable to study the effect of strong magnetic field and corresponding 
quantum corrections on compact neutron stars. 

There are strong reasons for believing that the hadrons are composed of quarks , and the idea 
of quark stars has already been existed for about twenty years. If the neutron matter density 
at the core of neutron stars exceeds a few times normal nuclear density a deconfining phase 
transition to quark matter may take place. As a consequence, a normal neutron star will be 



converted to a hybrid star with an infinite cluster of quark matter core and a crust of neutron 
matter. In 1984, Witten suggested that strange matter, e.g., quark matter with strangeness per 
baryon of order unity, may be the true ground state (Witten 1984). The properties of strange 
matter at zero pressure and zero temperature were subsequently examined and it was found that 
the strange matter can indeed be stable for a wide range of parameters in the strong interaction 
calculations (Farhi and Jaffe 1984). Therefore, at the core, the strange quarks will be produced 
through the weak decays of light quarks (u and d quarks) and ultimately a chemical equilibrium 
will be established among the participants. Since, the strange matter is energetically favorable 
over neutron matter, there is a possibility that whole star may be converted to a strange star. 

In the present paper, we study the effect of strong magnetic field on strange quark matter 
including interaction upto one order of strong coupling constant. Throughout the calculation, as 
a first approximation, we have assumed the existence of constant uniform magnetic field H 10 16 
G in the strange star. The interior of the quark star has high conductivity and has therefore 
zero electric field and uniform current density. The magnetic field H cx -j^ oc J ~ constant, 
where I is the total current, J is current density and R is the radius of the spherical quark 
star. Taking into account the low - temperature corrections in the thermodynamic potentials, we 
investigate the stability of strange matter for finite temperature, finite chemical potentials and 
finite magnetic field. At the end, we apply strange quark matter equation of state to structure 
equations of a relativistic spherical static star to calculate the structure parameters of strange 
quark stars. 

For a constant magnetic field along the z-axis (A = (Hy, 0, 0)), the single energy eigenvalue 
is given by (Landau and Lifshitz 1965) 



where n=0, 1, 2, being the principal quantum numbers for allowed Landau levels, s = ±1 
refers to spin up(+) and down(— ) andp, is the component of particle (species i) momentum along 




(1) 



the field direction. Setting 2n + s + 1 — 2v, where u — 0,1, 2..., we can rewrite the single particle 
energy eigenvalue in the following form 

Si = ^jvi 2 + W + 2uq t H . (2) 

Now, it is very easy to see that v — state is singly degenerate, whereas, all other states with 
v 7^ are doubly degenerate. Thus we set b u = 2 — 5 u q. For v — 0, b v — 1, which is the 
lowest Landau level. Since the temperature T << /Xj at the core of quark star, the presence of 
anti-particles can be ignored. Now instead of infinity the upper limit of v sum can be obtained 
from the following relation 

Pfi 2 = Hi 2 - m 2 - 2uqiH > 0, (3) 
where ppi is the Fermi momentum of the species i, which gives 

v < l — 1 = Vma^ {nearest integer). (4) 
Iqiti 

Therefore, the upper limit is not necessarily same for all the components. As is well known, 
the energy of a charged particle changes significantly in the quantum limit if the magnetic field 
strength is equal to or greater than some critical value = mi 2 c 3 /(qfi) (in G), where rrii and 
qi are respectively the mass and the absolute value of charge of particle i, h and c are the reduced 
Planck constant and velocity of light respectively, both of which along with Boltzman constant 
k B are taken to be unity in our choice of units. For an electron of mass 0.5 MeV, this critical 
field as mentioned above is ~ 4.4 x 10 13 G, whereas for a light quark of current mass 5 MeV, this 
particular value becomes ~ 4.4 x 10 15 G. 

Then the thermodynamic potential in presence of strong magnetic field H(> H^ c \ critical 
value discussed later) is given by 

A* = - 9 -^§^ J de-f: K^- ln[l + exp(^-e i )/T], (5) 

where gi is the degeneracy of the species i. 



Integrating by parts and substituting 



Pi = ± V 4 - m i - 1vq t H , (6) 

for all T, one finds 



__g_^H_ r ^ 2^e 2 -m 2 -2vq i H 

An 2 J ^ ^^ [exp{£t _ fIi)/T+1] V) 

where the sum over v is restricted by the condition e > y/m 2 + 2uqH and the factor 2 takes 
into account the freedom of taking either sign in eq([]). For T =0, therefore, approximate the 
Fermi distribution by a step function and interchange the order of the summation over v and 
integration over e, 

CV = bJ" de^ej - m? - 2 m H 

g iqi H 
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'raf + 2uqiH 
which is first order thermodynamic potential. 

Considering interaction upto one order of strong coupling constant a c (Freedman and McLer- 

ran 1978, Alcock et al. 1986 and Haensel et al. 1986), the second order thermodynamic potential 

in presence of strong magnetic field for zero temperature is given by 

2 g iqi H 8a c ^ J 2 2 /i, 2 - m t 2 - 2^ff | 

a =-^v?M ( ^~ m? * } ^ J ( ) 

The low - temperature corrections to the first order and second order thermodynamic potentials 
have been derived in Ref. (Sahu 1995) and these are: 
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This expansion is valid for , = <C 1. This means that the distance from the edge 

of any Landau level e v (p z ) = Jm% + 2qiHu to the Fermi surface fii is much greater than the 
temperature, or in otherwords, the Landau level with v = v max should be partially filled, i.e., its 
edge cannot concide with the Fermi surface. 

So, the total thermodynamic potential is 



&i(T, H, Hi) = lV(T, H, fii) + £V(T, H, fa). (12) 

In our study, we assume that strange quark matter is charge neutral and also chemical 
equilibrium, then 

A*d = lis = A* = y-u + He, (13) 
and charge neutrality conditions gives 

2n u — rid — n s — 3n e = 0. (14) 

The baryon number density of the system is given by 

n B = -(n u +n d + n s ). (15) 

Using above eqs(0, 0, |15|) one can solve numerically for the chemical potentials of all the 
flavors and electron. Having the expression for the total thermodynamic potential, we may move 
forward to calculate the number density (—dQi/d/ii) and the magnetization (—dQi/dH) of the 



species i (u, d, s, e) (Sahu 1995). Using Eq. dl2|) and definition of number density, one has the 
following expression: 
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Now, for T = 0, we have the number density of the species i: 
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In figure 1, we have plotted the number density of u quarks as a function of magnetic field for 
fixed chemical potential (fi u = 20 MeV). The u quark density is showing an oscillating behavior 
as consecutive Landau levels are passing the Fermi level. The number density is low for case (a), 
where a c = 0.0 and the number density increases with increase in a c = 0.01 for case (b) and 
a c = 0.05 for case (c) respectively. Thus, we noticed that the number density is high for strongly 
interacting quarks. 

In the weak magnetic field limit (Hq (/i 2 — mf)), one may reduce the number density as 
follows: 



n,(# «Gu 2 -m 2 ),/i,T = 0) 
~ JL {-(a 2 - m 2 ) 3/2 
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where, we can substitute in the limit 2q i Hi / = 9 



2 qi H 

E 

i/=0 



d v — lim 



1 



tf^o 2qiHJo 



(ii?— m?) 



d 



(19) 



The ratio of u quark density, n u (H <C (/i 2 — m 2 ))/n u (H) as function of chemical potentials for 
fixed magnetic field, curve (a) H xs 5 x 10 15 G and curve (b) H 10 16 G, is shown in figure 2. 
We choose the weak field limit, H to be 4.4 x 10 13 G. As the chemical potential increase, the 
ratio decreases and reaches to unity for both curves (a) and (b). We noticed that the magnetic 
field has significant contribution to the quark densities and hence to the total thermodynamic 
potential. For illustration purpose, we choose here the u quark density. 

on, 



From the definition of magnetization M = 
M(H) is the vacuum magnetization, one has 



dH 



, where M(H, fi, T) = M(H) + M(H, /i, T), 



M(H,n,T = 0) = 
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One may calculate the magnetic susceptibility x, which is defined as (dM/dH). 

Figure 3 shows the magnetization as a function of magnetic field for fixed chemical potential 
(fi u = 20 MeV). We have not considered the vacuum magnetization in this figure, because the 
vacuum contribution is small (Sahu 1995). We noticed that the quark gas exhibits the de Hass 
- van Alphen effect. We considered the strong coupling constant in curves (a) a c = 0.0, (b) 
a c = 0.01 and (c) a c = 0.05 respectively for illustration purpose. The magnetization is high at 



small value magnetic field and approaches to case (a) (e.g., free quark gas) for large magnetic 
field. Thus the de Hass - van Alphen effect amplifies with increase in strength of the strong 
coupling constant. 

The total energy density and the external pressure of the strange quark matter is given 
respectively by 

go. 

e = E^ + 5 + £r^-r(— f) w 

i i 

p = -J2^i~B, (21) 

i 

where i — u, d, s, e. The last term is the total energy density is entropy of the system, which 
is non zero and is given by 
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Here, we have considered the conventional bag model for the sake of simplicity in presence of 
magnetic field. We are assuming that the self interacting quarks are moving within the system 
and as usual the current masses of both u and d quarks are extremely small, e.g., 5 MeV each, 
whereas, for s-quark the current quark mass is to be taken 150 MeV. We choose the bag pressure 
B to be 56 MeV fm~ 3 and the strong coupling constant a c < 1. Also, we set the magnetic field 
to be H ~ 10 16 G in our calculations. Since we choose the magnetic field along the z-axis, it 
follow from energy-stress tensor that the constant uniform magnetic field (H = 10 16 G) contribute 
to the pressure and the energy density by |^-, which is much smaller than the bag pressure, e.g., 
0.0025 « 56MeV /m" 3 . 



The variation of energy per baryon with baryon number density is shown in figure 4. The 
curve (a) is for free quark gas, a c = 0.0 at zero temperature and curve (b) is for interacting quark 
gas with interaction strength a c — 0.1 at zero temperature, whereas, the curve (c) is for a c — 0.1 
and temperature T = 50 MeV. For all the cases, the magnetic field strength is 10 16 G. We 
concluded from the figure that the interaction strength and finite temperature make the strange 
quark matter energetically more stable compare to the zero temperature and zero interaction 
strength. 

Next, we have shown the variation of pressure with energy density in figure 5. These are the 
equation of states of strange quark matter. The curve (a) is for free quark gas (a c = 0.) at zero 
temperature and curve (b) is for interacting quark gas (a c = 0.1) at zero temperature. From this 
figure, we noticed that the interaction strength makes the quark equation of state soft. Also, we 
have seen that the softness increases with increase in interaction strength. 

The mass and radius for nonrotating strange quark stars are obtained by integrating the 
structure equations of a relativistic spherical static star composed of a perfect fluid which is 
derived from Einstein equation. These equations are given in Ref. (Datta et al. 1992, Sahu, 
Basu and Datta 1993 and Ghosh and Sahu 1993), hence, we are not reproducing here. For a given 
equation of state, and given central density, the structure equations are integrated numerically 
with the boundary conditions m(r = 0) = 0, to give R and M. The radius R is defined 
by the point where p ~ 0. The total gravitational mass M, moment of inertia I, surface 
red shift z and the period Pq corresponding to fundamental frequency Qq ar e then given by 
M = m(R), I = I(R), z = (1-2GM/Rc 2 )~^ 2 -1 andP = g respectively, where fi = (^) 1/2 
(Cutler, Lindblom and Splinter 1990). These are correspond to maximum mass of stable star 
that can have, presented in Table 1. Figure 6 shows the variation of mass with central density 
for two equation of states as illustrated in Fig. 5. We noticed from this figure and table that 
with increase in interaction strength, the maximum mass and the corresponding radius decrease 



from 1.83 to 1.64 solar mass and from 10.4 km to 9.58 km and therefore, the star becomes 
more compact. Similarly, the values for surface red shift, moment of inertia and fundamental 
period decrease. Also, we have noticed that with increase in magnetic field strength, the density 
also increases and therefore, the star becomes more compact (Chakrabarty and Goyal 1994, 
Chakrabarty and Sahu 1995 and Chakrabarty 1995). 

In conclusion, we concluded that the presence of strong magnetic field in interacting strange 
quark matter reduces the mass and radius of strange star. That is the star becomes more compact. 
With inclusion of interaction strength of quarks at finite temperature and at zero temperature in 
presence of strong magnetic field, the strange quark matter is more stable compare to free quark 
gas at zero temperature. 
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Table 1: The radius (R), mass (M), surface red shift (z), moment of inertia (/) and period of 
fundamental frequency (P ) of strange stars versus central density e c for two cases; (a): a c = 0.0 
and (b): a c = 0.1. For both the cases H m 10 16 G. These values are correspond to maximum 
mass of stable strange star. 
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Figure 1: The number density of u quarks as a function of magnetic field for fixed chemical 
potential (//„ = 20 MeV). Where Case (a) is for a c = 0., case (b) is for a c = 0.01 and case (c) is 
for a c = 0.05 respectively. 

Figure 2: The ratio of u quark density, n u (H <C (fj% — m^))/n u (H) as function of chemical 
potentials for fixed magnetic field, curve (a) H ks 5 x 10 15 G and curve (b) H ks 10 16 G. The 
choice of the weak field limit, H is 4.4 x 10 13 G. 

Figure 3: The magnetization as a function of magnetic field for fixed chemical potential (fj, u = 20 
MeV). The curves (a) is for a c = 0.0, (b) is for a c = 0.01 and (c)a c = 0.05 respectively. 

Figure 4: The variation of energy per baryon with baryon number density with constant 
magnetic field H ks 10 16 . The curve (a) is for free quark gas, a c = 0.0 at zero temperature and 
curve (b) is for interacting quark gas (a c = 0.1) at zero temperature and curve (c) is for a c = 0.1 
and temperature T = 50 MeV. 

Figure 5: Pressure and energy density (e) curves for case (a) a c = 0.0 and case (b) a c = 0.1 
with constant magnetic field H m 10 16 . 



Figure 6: Mass and central density (e c ) curves for two cases as mentioned in figure 5. 



